March 4, 2013 12:49 WSPC/INSTRUCTION FILE Article-final 



Matrix Product State representation for Slater Determinants 
and Configuration Interaction States 



o 

PIETRO SILVI 

^ Institut fur Quanteninformationsverarbeitung, Universitdt Vim, D-89069 Vim, Germany 

International School for Advanced Studies (SISSA), Via Bonomea 265, 1-34136 Trieste, Italy 



DAVIDE ROSSINI 

__ NEST, Scuola Normale Superiore and Istituto di Nanoscienze - CNR, Pisa, Italy 

Oh 



ROSARIO FAZIO 

NEST, Scuola Normale Superiore and Istituto di Nanoscienze - CNR, Pisa, Italy 
Center for Quantum Technologies, National University of Singapore, Republic of Singapore 

GIUSEPPE E. SANTORO 

International School for Advanced Studies (SISSA), Via Bonomea 265, 1-34136 Trieste, Italy 
CN International Centre for Theoretical Physics (ICTP), P.O. Box 586, 1-34014 Trieste, Italy 

^ CNR-IOM Democritos National Simulation Center, Via Bonomea 265, T34136 Trieste, Italy 

in 

y—i VITTORIO GIOVANNETTI 

NEST, Scuola Normale Superiore and Istituto di Nanoscienze - CNR, Pisa, Italy 

in 
o 

£N| Slater determinants are product states of filled quantum fermionic orbitals. When 

T~H they are expressed in a configuration space basis chosen a priori, their entanglement is 

!* bound and controlled. This suggests that an exact representation of Slater determinants 

. as finitely-correlated states is possible. In this paper we analyze this issue and provide an 

exact Matrix Product representation for Slater determinant states. We also argue possible 
j_j meaningful extensions that embed more complex configuration interaction states into the 

description. 
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1. Introduction 

Starting from its first formulation in 1992, the density matrix renormalization group 
(DMRG) has rapidly established itself as the leading numerical method in the sim- 
ulation of statical and dynamical properties of strongly correlated one-dimensional 
quantum systemsSH The versatility of this approach allowed the extension of 
DMRG applications to fields quite far from those it was originally designed to 
address, e.g. the quantum chemistry of small to medium-sized molecules^. A deeper 
understanding of the mechanisms which contribute in making DMRG such a success- 
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ful tool has been achieved recently, when this technique was shown to be equivalent 
to a variational approach over the set of Matrix Product States (MPSsjl The lat- 
ter were originally named "finitely correlated states" , due to their intrinsic ability 
in capturing the internal structure of many-body quantum states having exponen- 
tially decaying correlations^. Remarkably it has been found that the ground state 
of the one-dimensional AKLT model is exactly given by a MPSS1Z| thus stimulating 
a plethora of investigations in the quantum many-body realm (e.g. see Refs. [H [9l 
ITUl [TTJ [12j [13l HU [15] and references therein) . The MPS architecture is so imme- 
diate and powerful that in some cases it is even amenable to analytical studies, 
providing an exact, and often optimal representation for several classes of quantum 
states beyond the AKLT model, including GHZ, W, Majumdar-Gosh, Cluster^!, 
and Laughlin states^. 

In this paper we discuss how Slater determinants can be represented in terms 
of MPSs. Slater determinants correspond to product states of (say) N fermions 
where each fermion occupies a different single-particle orbital out of an orthonor- 
mal set (typically solution of a mean-field Hamiltonian) . They are the starting point 
of most many-body calculations (e.g. ab-initio Hartree-Fock and post-Hartree-Fock 
methods)?^ and, even though lacking true electron-electron correlations, their char- 
acterization poses significant problems when working on a configuration space for 
which the single-particle orbitals entering the determinant arc not local (by this we 
mean that the latter might have extended and non-mutually excluding supports in 
a wavefunction basis chosen a priori). For this reason, most previous attempts to 
address quantum chemistry problems with DMRdHEH, and more in general with 
Tensor Network variational structures' 20 ! 21 !, agreed on adopting the Hartree-Fock 
basis as the selected configuration space for the correlation calculus. In the present 
work, we step back from this assumption, and show that that the description can be 
kept simple and manageable even when the chosen configuration basis is not related 
to the solutions of the one-body problem. 

The paper is organized as follows: in Sec. [2] we review the basic concepts of 
fermionic quantum states and set the notations. In Sec. [3] we provide our prescrip- 
tion to describe a generic Slater determinant state by means of a MPS. To this 
aim we first introduce a representation of a single fermionic operator in terms of 
a Matrix Product Operator (MPO) and then exploit that design in order to build 



the Slater determinant MPS (Subsec. 3.2). In this framework it will then be easy to 



introduce a compact and numerically manageable representation for the many-body 



transformation which corresponds to the one-body basis change (Subsec. 3.3). In 
Sec. [4] we will prove the optimality of our description in a general setting, through 
some considerations based upon quantum entanglement estimators. A straightfor- 
ward generalization to other particle-exchange statistics is given in Sec. |5j while 
a possible extension to Configuration Interaction is described in Sec. [6] Finally in 
Sec. [3 we draw our conclusions. 
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2. Basic notations 

Consider a system composed of N spinless fermions, which can occupy L states 
(sites), identified by a canonical basis of one-body wavefunctions, in some effec- 
tive configuration (or momentum) spaced As usual we can account for the Fermi 
statistics of the system by defining a complete ordering of the canonical basis and en- 
forcing a Jordan Wigner Transformation (JWTf W 4 \' 2b lB, which is common ground 
when addressing fermionic problem with numerical renormalization strategies^. 
Accordingly, each lattice site, say the one labeled with the index £, is fully charac- 
terized by assigning a set of Pauli matrices a" (a — x,y,z), of = (af ± icr|)/2, 
and indicating with |0)^, |1)^ its empty /occupied level, which satisfy the relations 
cr| \0) e = \0) e , <j\ \X) t — — |1) g respectively. Then, following the standard procedure, 
we construct the correspondence: 

|n>— Ho> 1 |o) a ...|o) £ = |oo...o) , (i) 

ct — >■ eg = of ® . . . ® o-JLj ® of ® ® . . . ® t L , 
c\ct — > c}ct = (1 - of)/ 2 , 

where |0) is the vacuum (Fock) state of the fermionic system, while c\ (q) is the 
creation (destruction) operator which creates (annihilates) a fermion on the £-th 
one-body state of the canonical basis, and obeys the canonical anti-commutation 
relations {cj, cj„} = {eg, c m } — 0, {q, cj^} = <^, m . As long as the set of L one-body 
levels is complete, any Fermionic state \^} can now be expressed as a many-body 
state f l qubits via the following mapping: 

i*)= E ■■■ E r sl ... s A4r.-.(clr L P) 

si£{o,i} s L e{oa} 

_^ l$(^)= ^ Tsi ... S£ | Sl ... Si ), (2) 

{s 3 }e{o,i} 

where, in the first line, (ct)° stands for the identity operator 1, while the creation 
operators c\ are placed according to the chosen ordering. The tensor T Sl ... SL holds 
the quantum amplitudes of every fermionic occupancy configuration. 

A Slater determinant state |E) is a special iV-particle state that can be cast in 
the form 

\z) = c{%..J N \n), (3) 

where the are creation operators associated to a collection of mutually orthogonal 
one-body orbitals. These are eventually delocalized with respect to the original 
canonical basis. In particular, we decompose the c^s according to 

L 

4=5>«(*)ct> (4) 
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Figure 1. (Color online). MPS representation of the tensor Ts 1 ...s L appearing in Eq. pj. Every 
element A^' is a three-index complex tensor, graphically depicted as a red box. The two yellow 
elements represent the correlation boundary vectors. Open links, the vertical ones, refer to physical 
degrees of freedom (sites); while the closed, horizontal links, are the correlation space indices, which 
are contracted through the matrix product. 



with 4> a (£) satisfying the orthonormality relations 



L 



£=1 



(5) 



The explicit determinant form of the many-particle state |E) becomes manifest when 
expanding the product c\c\ . . . cjy in terms of the cjs. Indeed simple mathematical 
manipulations allow one to write 



IS) 



E 

1<<i 



<L 



M*i) ■■■ M*n) 



(6) 



which can be cast in the form of Eq. ([2| by taking 



E 



i<e 1 <...<e N <L 



n ^ 

>5 



(7) 



where t' spans over the site indices which appear in the determinant, and £" on those 
not appearing. Equation ([7]) explicitly manifests the complexity of representing | E) 
in qubit-oriented configuration space. 

The purpose of the paper is to provide an exact and efficient representation for 
states of the form ([6]), |7| in the MPS formalism. In practice this consists in writing 
T Sl ... SL as a contraction of smaller elements according to the decomposition 

T Sl ... SL = (b \AW-AW-...-AM\b L ), (8) 

where As are a set of D x D matrices, while and (6q| are respectively a D- 
dimensioned column vector, and a /^-dimensioned linear functional (row vector JP^. 
A pictorial representation of this identity, which follows the standard graphical 
convention^, is presented in Fig. [l] In Eq. ^ D is a free parameter quantifying the 
degree of refinement of the MPS representation (thus also referred to as refinement 
parameter, as well as bond dimension). As detailed in Sec. [4j we will focus on an 
explicit construction ^ which ensures the smallest possible value of D, regardless 
of the shape of the orbitals <f> a (£) entering in Eq. (JT]). We will show that this minimal 
(or optimal) bond dimension is D = 2 N . 



March 4, 2013 12:49 WSPC/INSTRUCTION FILE Article-final 



MPS representation for Slater determinants and CI states 5 



3. MPS representation of a Slater determinant 
3.1. MPOs for delocalized Fermi operators 

To find a MPS representation Q for |S) it is instrumental to first provide a MPO 
description^ for the Fermi operators introduced in Eq. Q. According to this goal, 
for each lattice site I and for each orbital a, wc identify a collection of four D x D 
matrices B^ , B^ , B^\,B^\ and boundary D-dimensional vectors (/3 | and |/3l), 
such that the creation operators can be decomposed as 



1 

E 

Sl...«t=0 r!...r L =0 



E 



mzKr ■■■■■ B [L] ? L \h) (c\r ■ • • (4n wi # 



The tensor 



■ ,sl 



0. 



(9) 
(10) 



entering the above expression is the MPO representation of depicted in Fig.|J At 
the level of the JW mapping ^ this correspond to decompose the anticommuting 
qubit-orbital construction operators in the following formalism: 



HP 



E E (fti* 

»i...ai,=0 ri...rz,=0 
where such operators are defined by 

4 = E^)4 = E^w 



.B[ i ^| / 8 £ ),|« 1 ... ai )<r 1 ...r z | 



(8 1; 



(11) 



(12) 



Elaborating from the results of Ref. [30l we claim that to build such representation 
it is sufficient to take D = 2. Indeed, assume to add a two-level fictitious degree 
of freedom, with canonical basis states |0) and |1). Then it is possible to write the 



previous expression ( 12 1 as a contraction over this new degree of freedom of a matrix 
product: 



(II 



— }Li 

[] ]±t ® |0)(0| + of <g> |1)(1| + M£)aJ <g> |1)(0|] |0) 



(ii n 

t=li- 

(13) 







Bp] 




fl[3] 


• • • 







Figure 2. (Color online). MPO representation of the fermionic operator cj,, as appearing in 
Eq. {9}. Every blue box i?^' is a four-index complex tensor. 
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where the product is sorted in increasing order in I from left to right. The Be terms 



in Eq.(13) are meant as matrices in both the £-th site and the common fictitious 



degrees of freedom, and read 
1 



B e = - (l t + o\ ) ® (1/ + of) 



(14) 



To recover the formalism in (10), we just need to write Bi entrywise in the real 
degree of freedom. When doing this we obtain 

bK = t f 



Bl% = o 



) = (J)=I°) O8ol = (0 1) = (11- 



(15) 



It is easy to see that the information on the orbital a, via the orbital wave- function 
4> a {£), enters only in one of the 16 elements of the four-index tensor Bf. apart from 



that, the expression (15) is formally homogeneous in i and no dependence on a is 
left on the boundary elements (/Jo| and |/3l). 

In a similar way we can also construct the MPO representation for the annihi- 
lation operators Cf.. Alternatively we find it useful to recall that the MPO represen- 
tation behaves very well under the adjoint application: in practice, if the matrices 
B^ s provides a MPO representation for the operator O then B*^ r do the same for 
the adjoint operator 0t (where x* denotes complex conjugation of x), i.e. 



l 

E 

s 1 ...s L =Q r!...r L =0 



E 



(MB* 



B 



\p L ) (4r...(4nfi><oi 



,t \»L\ 



As an final remark, we notice that, while still keeping the matrices B^' 



(16) 
enter- 
ing ([10]) as in Eq. (Il5j) it is possible to modify the action of the associated operator 
by simply changing the boundary vectors (/?o| and (0l\- For instance, by identifying 
both | /Jo) and with |0), the resulting MPO will represent the identity operator. 
More generally, manipulating boundary vectors gives 



(A>l 


\h) 


operator (spin) (Fermi) 


(name) 


(11 


|0) 


ft 

Ol L a 


Construction 


(0| 


|0) 




Identity 


dl 


11) 


(o z )® L (l-2Jcf 


,L Parity 


(o| 


11) 





Null. 



(17) 



In particular, when = |0) one can then regard the left correlation space bound- 
ary as a local switch that "activates" or "deactivates" the creation operator. 



3.2. Stacking Fermi operators into a MPS 

Let us now focus on the Slater determinant state |S) in Eq. following a prescrip- 
tion analogous to the one introduced for MPS representations of algebraic Bcthe 
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AM 



r [ 7 i 



9 



i — — — rj— 



- .H 



T 



9 



Figure 3. (Color online). The MPS construction provided in Eq. | |19| is obtained by vertically 
stacking Fermi operators in their MPO form. Thus an tensor is the contracted column of 
tensors applied to the local vacuum state |0). 

ansatJ^. The vacuum state |fi) — > |0...0) is a MPS with D = 1, since it is a 
product state. It is clear that, when we apply to a given MPS |$mps)i we can 
exploit its MPO formalism presented in the previous subsection, and immediately 
obtain a MPS expression for |$') = ctJ$ MPS ). A sketch of the procedure is depicted 
in Fig. [3] 

This process can be repeated recursively, every time adding a new & a until the 
whole operators string c\ . ■ ■ cjy is included, and acting on |fi). This is equivalent 
to writing the following three-index tensor, which represents the stacking of B^ 
tensors applied to the local vacuum coming from the Fock state 



92— 9w=0 



' B [41]« 8B [i.2]*3.. 



9BV^ J) , (18) 



here the B^'°^ r s compose the set of 2 x 2 matrices introduced in Eqs. ( 15 1, concerning 
the orbital cj> a . This description uses as a whole a total refinement parameter (also 
called bond- link dimension) of D = 2 N ; exponential in the fermion number, but 



disregarding the total size L. By substituting (15) into Eq. (18), and exploiting the 



fact that Bq is always the null matrix, the whole equation is reduced to 



A = 1 2 n x2 n 



N 



(19) 



The correlation boundary vectors are similarly constructed: 

\b L ) = \P L r N and (6o| = (A)| 8W , (20) 
i.e. \b L ) = |0)® w = (10 0---0) T and (b Q \ = (1\® N = (0---001). With the pre- 



scriptions introduced in Eqs. (1 1 91) and (20) and a little algebraic manipulation, the 
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reader can check that such MPS provides the correct amplitudes for the Slater 
Determinant state 

l 

|E> = (b \A^-...-A^\b L )(c\r...(c[r-\n). (21) 

Sl...SL=0 



Example - We would like to end this section by showing a simple example 
about this representation at work. Let us consider the case of N = 2 fermions 
occupying two distinct orbitals. In this scenario A$ is the 4x4 identity operator, 
while Af 1 = <t>i{t)[<J- ® 1] + 4>2{£)[<J Z ® cr-], i.e. 

/ 0\ 



N = 2 



Af = 



(22) 



W) o oo 

o 

V o 0r {i) -M*)o/ 

The only products of matrices leading to nonzero amplitudes are those where exactly 



two excitations |1) are present. Thus the sum in Eq. (21 1 reduces to 



Mh) 0i (£2) 



E 



(23) 



where we have explicitly recovered the determinant expression, thus confirming the 



validity of the MPS construction ( 19 ) and ( 20 ) 



3.3. 



Tensor network representation of one-body wavef unction basis 
change 



When we derived the MPO representation (151 for & a , we also mentioned that it is 
possible to control its overall action by adjusting one of the two correlation boundary 
vectors, say the left one (/3q| [see Eqs. (17)]. In particular, the MPO produces b\ 
if (/3o I — (1|, while it coincides with the identity operator 1 as we set (/3o| = (0|. 
Recall that the fermionic orbitals Q (£) we adopted for the Slater determinant state 
formed an orthonormal set: let us complete it to an orthonormal basis {<fi a (£)}a, 
with a £ {1 . . . L}. The dimension must be L by the assumption that the original 



set of L wavefunctions was complete. As before, Eq. ( 15 ) provides the construction 
operator cj, for any of those orbitals (f) a (£). 

Now we stack together the MPOs, like we did for the Slater state, but instead 
of using only N of them, we stack the complete set, ordered from a = 1 (on top) 
to a = L (at the bottom). Moreover, instead of using the standard left correlation 
boundary vector (/Sola = (1| we se t a generic (/3oU = {q a \, where q a are classical 
binary variables. It is clear that the many-body operator O generated from this 
setup is equivalent to 



Q=~c\ qi c\ q2 ct 9L . 



(24) 
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Figure 4. (Color online). Tensor network representation for the many-body state transformation 
of one-body basis change. 7~ old and 7~ ncw from Eq. ( |26| are tied through this pictorial equation. 
Blue boxes correspond to B^ l '°^ tensors, where I and a respectively label the horizontal and vertical 
axis position in the grid, starting from the left-bottom corner. 

Finally, we apply such operator to the vacuum |fi). By doing so we are actually 
defining an application from the binary strings of {q a } a to the Fermi system: 

{ qi ... qL \^\^t)=A qi tl q2 ---A qL \n)- (25) 

By linearity, this map extends to all the space generated by [q\. . .qi\, which cor- 
responds to the whole correlation bond-link space, {qi...qi\ being its canonical 
product basis. The map is clearly invertible, and its inverse is basically a JWT from 
the Fermi space to the qubit lattice. Only, this time the extended orbital set <fi a has 
been chosen as basis of one-body wavef unctions, so it is formally similar to ([!]), but 
the one-body basis has changed (the old one is associated to the c, the new one to 
the c). 

In conclusion, we could use the extended MPO stack formalism to represent a 
transformation of the many body state \^f) corresponding to a change of the chosen 
basis of one-body wave functions. I.e. assuming that we can write in both of the 
following expansions 

l l 

i*>= E v! d . SL (c\r...(cir-\n)= £ r^^r-.-^n^ (26) 

si...s_l=0 q\... gi=0 

then the two components tensors T old and 7~ ncw satisfy the relation depicted in 
Fig-i 

As a concluding remark, we recall that 2D tensor lattices cannot be efficiently 
contracted exactly in general. Nevertheless, the tensor network that we built in 
Fig. Q is made of sparse tensors (the high amount of zero tensor elements due 
to particle-conservation symmetry relations), thus effectively making it an easily 
contractible network. 
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4. Entanglement and optimality of the description 

The mere existence of an exact finitely-correlated-state representation for a generic 
Slater Determinant |£), automatically defines an upper bound to the entanglement 
contained in |S). Indeed, MPS manifest a well-known bounds™ to their entanglement 
entropy (the von Neumann entropy of a subsystem-reduced density matrix p) given 

by 

Svn Go) = -Tt\p log p] < log D. (27) 

Clearly, any class of states allowing an exact MPS representation with established 
(eventually non-scaling) bond dimension D, undergoes the same upper bound. 
We now argue that, if no further information upon the filled orbitals <j) a {£) within 



|£) is assumed, the exact representation we just gave, Eqs. ( 19) and (20 1, is the most 
efficient in terms of MPS. By this, we mean that there cannot be another general 
MPS representation A for arbitrary |S) spending a bond- link dimension D < 2 . 
To prove this, we will consider a specifically built ./V-fermions Slater Determinant, 
and show that its entanglement entropy is equal to <Svn = N log 2 (here we are 
assuming that L > 2N); since any exact MPS representation for this state should 
have D > 2 N due to (27), we must conclude that A fails to reproduce this state, and 



hence it is not general. The present entanglement bound estimator can be perceived 
as a special case of those treated in Refs. l52"HT7l but it was developed autonomously. 

For this setting, let us adopt a special set of (L/2-periodic) plane- waves 4>a(£) = 
L~ 1 / 2 exp(4iritt£/L). Now we consider a half-system partition, and define a new 
double set of orbitals {(ffflil), <p^(£)} a respectively reducing the original cf) a (£) to 
the left half and right half supports: 

= V2 Q(L/2 - £) cj> a (£) and <f>W(£) = V2 Q(£ - L/2) <j> a {£), (28) 

with O being the Heaviside step function. Even though in a general case a new set 



of wavefunctions generated via ( 28 ) would no longer be orthonormal, this specific 



choice of original orbitals <p a (£) preserves orthonormality: 

1 (29) 

E ^ {£) = E € ] wf (£) = s a , . 

Let us write Fermi operators corresponding to this new set, satisfying the anti- 
commutation rules {c a .L, Cg L } — {c a ,R, Cg R } = S a ^, all the other anticommu- 
tators being null. It is clear that the original c decompose in the new ones as 

thus allowing us to write the whole Slater determinant 

state as: 



l S ) = ^72 (A,L + 4, R ) ■ ■ ■ (4,L + 4,il) (30) 
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By partially tracing over degrees of freedom related to one half of the system, say the 
right one, we achieve the reduced density matrix of the left half p\ = Ti R [|£)(£|]. 
With this goal in mind, just set |S') = ci|£) and consider: 



1 



pi = - Tr fl [(£{ L + cl R ) |E')<E'| {c 1>L + c hR ) 

= \ (i* iL Tr R [|E')<E'|] c 1<L + Tr R [|S')(S'|]) = - (c{ L pf c h L + (% 



(31) 



where we exploited the cyclicity of the trace over right support operators c atR , and 
the fact that ci^JE') = 0. Also it turns that p\ and (c[ L p\ Ci t i) have orthogonal 
supports and the same spectrum^. Therefore we have iSvn(Pl) = Svn(Pl ) +log2 
[see pag.513 of Ref.H Theorem 11.8 (4)]. 



Finally, we repeat the argument (31) on |E'), and proceed by induction. In 



conclusion, we can claim that pfj is (isometrically equivalent to) 2~ N t 2 N x2 N i the 
maximally mixed state on a 2 N dimensioned space, having a von Neumann entropy 
of N log 2. This concludes the proof. 

The previous result allows a clear and sensible interpretation: fermions occupy- 
ing the various orbitals must be mutually uncorrelated due to the Slater determinant 
state nature, so the available entanglement is given by the self-correlation of every 
orbital, separately accounted. Rephrasing, there is no many-body entanglement in 
|E), but one-body entanglement of each particle can still be present. A single delo- 
calizcd fcrmion can contribute to the total amount with at most the entanglement 
of a unit (i.e. the amount of entanglement shared by a spin singlet), so N units is 
naturally the overall maximum. 

5. Other exchange statistics 



It is important to remark that the a z matrix in Eqs. ( 15 ) is the only responsible for 



establishing the correct anticommutation relations of Fermi statistics. That said, it is 



arguably easy to adjust (15) in order to adapt the representation to non-fermionic 
statistics^, as long as the particles involved are hard-core. When a single level 
cannot be occupied by more than one particle (guaranteed by the Pauli exclusion 
principle for fermions, must be inquired through other means otherwise) then it 
is again possible to map the physical model into a ID qubit lattice: the JWT of 
Eq. ([T]) actually becomes 

|0) — ► |00...0) (32) 

Q — > Wl ® . . . <g) Wf_i <g) (Tj <g> 1« <g) . . . ® t L , 

where W is the particle-exchange matrix which encodes the statistics. The 2x2 
matrix W is equal to the identity 1 for bosons, to the Pauli matrix a z for fermions. 
For abelian anyons which acquire a fractional phase tp under exchange, W is still 
a diagonal matrix^!, precisely it becomes the phase gate W = exp(«</j<7 2 ). Finally 
W can assume any shape when we consider particles whose exchange leads to the 
application of a non-trivial operator, that is non-abelian anyons^. 
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The generalization of the MPO representation (151 of the creation operators is 
quite straightforward for bosonic and abelian anyonic scenarios. One needs only to 
change the matrix as 

B [e] \=W (33) 

and keep the rest. Since now W is the phase gate (with trivial phase ip = if we are 
considering bosons), it is still diagonal in the eigenbasis of a z , thus meaning that 
the matrix product contraction with the same boundaries (/3q| = (1| and = |0) 



is reduced like in Eq. (12 1, i.c 



<t>a(Z) [VWi ® . . . <g> Wi-i ® a' ® l e+1 ® U] — > Y <4 = 4> ( 34 ) 

with a] (a^) being the abelian anyonic/bosonic construction (destruction) operator 
on site £, while a) a (a a ) the construction (destruction) operator for the particle 
orbital a. We can then write the MPS representation A for the hard-core many- 
body state free of correlations (generalization of the Slater Determinant for an 
arbitrary abelian exchange statistics W = e %H>cr ) as follows 

A^q = l 2 iv X 2N (35) 



N 



Unfortunately, the non-abelian anyonic case is more difficult to be addressed: 
since the general form of a particle-exchange matrix W does not commute with a z 
the MPO construction cannot eliminate rigorously multi-particle creation events, 
therefore the previous generalization is not exploitable. In this case, it is not clear 
to us whether a larger refinement parameter D is required. 



6. From Slater determinants to configuration interaction picture 

More accurate results can be achieved by adopting the Configuration Interaction 
(CI) picture. In this description, Hartree-Fock solutions are adopted as a canonical 
vector basis of orbitals for further calculations, embedding entanglement between 
particles^. According to such approach, one is interested to express quantum cor- 
relations by superposing few to several Slater determinant states, which typically 
share most of the HF orbitals while differing for a limited amount of those, in order 
to keep the calculation manageable for practical purposes. 

Having this scheme in mind, we would like to extend our previous Slater MPS 
(via MPO stack) representation to embed also CI states, where different orbital 
excitations are coherently overlapped. The ultimate ingredient of this perspective 
would be writing exact Matrix Product representations for every operator generated 
by the Fermi ones through sums and multiplications. Unfortunately, providing 
a general, elegant and high-compatibility formulation is no easy task, since the 
matrix product representations are likely to be operator-dependant and definitely 



March 4, 2013 12:49 WSPC/INSTRUCTION FILE Article-final 



MPS representation for Slater determinants and CI states 13 

not unique. Here we will consider the simplest nontrivial case, and speculate on 
various proposals. 

For instance, consider the expression 



02+2 = a c\c\ 



f3c, 



3 L 4- 



(36) 



We are now going to describe 02+2 as a MPO. Depending on whether we focus 
on the adaptability of the description or the economy on the refinement parameter 
dimension we can end up with different representations. 

In the first case one exploits the standard technique to sum coherently Matrix 
Product objects, and is strongly based on the MPO structure of Fermi operators 



in Eqs. (15 1. It is highly suitable for further generalizations, but at the cost of a 



sub-optimal bond-link dimension. Let us now adopt D = 8 and consider 



£[£,2+2]* = 











B [<,3]* 8jB [M]J 



(37) 



where the B^'°^ tensors are those defined in ( 15 ) for cj,. The basic idea behind this 
construction is to use a correlation space which is the Cartesian sum of the two 
original correlation spaces, and a matrix product object which is the block diagonal 
composition. Similarly we define the correlation boundary vectors, which contain 
information on a and /3: 



o 




p 




w 



a 



f 1 } 




f 1 } 






















0/3 





-(z)° 





w 




W 







(38) 



where we used distributivity of the tensor product <£> with respect to the Cartesian 
sum ©. Similarly, (b \ = (n) (8 (o o o i). 

On the contrary one might also like to keep the lowest possible value of the 
refinement parameter, D = 6. In this case one shall take 



B [e,2+2] 



1 



6x6 ■ 



( o 

-afait) 

pcj> 4 (£) 
a<j) 2 (£) 

V o M 







o o\ 
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and 



,2+2] 1 _ 
1 — 



/l o 
-1 




\Q 






-1 
-1 



o o\ 







0-10 
1 



(39) 



while boundaries are left as before \bi,) = |0) and (&o| = (o ...o i) = (5|. By multiply- 
ing the £?^> 2 + 2 l matrices it is easy to see that we are reproducing the correct action 
of the operator, i.e. 

{a(<i>i(ii)Mt*)-Mh)Mtij) +p(Mti)Mh)-Mh)Mtij)} C 14 2 - 

(40) 

One can ask whether this Matrix Product representation is optimal in terms of 
correlation bond- link dimension. A state of the form 62+2!^) manifests an entan- 
glement entropy at most equal to 5/2. This implies that a faithful MPS description 
would require a D > \/32, so that D = 6 is the smallest allowed integer, and thus 
arguably optimal. 

The present proposal presents various options for generalization, although find- 
ing the analytical MPO expression for a generic operator which is cheapest in terms 
of D is definitely non-trivial. 



7. Conclusions 

In this paper we studied matrix product representations for uncorrelated many-body 
states, expressed in an arbitrary configuration space upon which the orbitals might 
be highly delocalized. To achieve this goal, we provided the exact representation of a 
(delocalized orbital) creation operator as an MPO, which we checked to be optimal 
in terms of the refinement parameter. Stacking these MPO leads to the description 
of the uncorrelated many-body states quite automatically. The representation we 
introduced is tailored upon fermionic exchange statistics, but it generalizes easily 
to particles obeying any abelian exchange statistic. Also, it leads to a compact and 
simple tensor network representation of the many body transformation correspond- 
ing to one-body basis change. Finally, we discussed some possible extension of the 
previous description to embed also states carrying a controlled amount of many- 
body correlation, such as CI states. We provided some proposals for the simplest 
cases and argued that there is no unique, highly-compatible and optimal solution. 

We conclude by arguing that, while the tensor network representation proposal 
we tailored in this paper was based on a one-dimensional qubit lattice, we are 
positive that similar construction can be adapted to other lattice geometries, as long 
as some peculiar topological features are preserved. In particular, we are thinking 
about finitely correlated states on Cayley trees^, sharing with ID OBC systems 
the fact that no closed loops in the geometry are present. 
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